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A RIAF Interpretation for the Past Higher Activity of the Galactic 
Center Black Hole and the 511 keV Annihilation Emission 
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Abstract 

There are several lines of evidence that the super-massive black hole at the Galactic center had higher 
activities in the past than directly observed at present. Here I show that these lines of evidence can 
quantitatively and consistently be explained if the mean accretion rate during the past ~ 10^ yrs has been 
~ 10^~^ times higher than the current rate, by the picture of radiatively inefficient accretion flow (RIAF) 
and associated outflow that has been successfully applied to Sgr A*. I argue that this increased rate 
and its duration are theoretically reasonable in the Galactic center environment, while the accretion rate 
suddenly dropped about 300 years ago most likely because of the shell passage of the supernova remnant 
Sgr A East. The chance probability of witnessing Sgr A* in such a low state is not extremely small (~ 
0.5%). The outflow energetics is sufficient to keep the hot (^ 8 keV) diffuse gas observed in the Galactic 
center region. Then, I show that a significant amount of positrons should have been created around the 
event horizon during the higher activity phase, and injected into interstellar medium by the outfiow. The 
predicted positron production rate and propagation distance are close to those required to explain the 
observed 511 kcV annihilation line emission from the Galactic bulge, giving a natural explanation for the 
large bulge-to-disk ratio of the emission. The expected injection energy into interstellar medium is ~ MeV, 
which is also favorable as an explanation of the 511 keV line emission. 
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1. Introduction 

A variety of active and high energy phenomena are seen 
in the direction towards the Galactic center over a broad 
range of wavelengths. It is well established that the cen- 
ter of our Galaxy, Sagittarius (Sgr) A*, is a supermas- 
sive black hole (SMBH) with a mass of M, ~ 3 x lO^M© 
(Genzel et al. 1997; Schodel et al. 2003), and a number 
of models have been proposed to explain the radiation 
and its spectral energy distribution (SED) from Sgr A* 
by accretion of matter onto it (sec Baganoff et al. 2003 
for a review and comparison to X-ray observations). The 
concept of radiatively inefficient accretion flow (RIAF) is 
theoretically well-motivated based on the physics of ac- 
cretion flows, and models based on the RIAF picture have 
successfully been applied to many low accretion rate sys- 
tems with M ^ O.IA/Edd, where A/Edd is the Eddington 
mass accretion rate (Narayan & Quataert 2005 for a re- 
view). Sgr A* is one of the best studied examples, and 
the model of Yuan, Quataert, & Narayan (2003, 2004, 
hereafter YQN03 and YQN04, respectively) can repro- 
duce the observed SED of Sgr A* in a wide range of wave- 
lengths from radio to X-ray bands. An important ingredi- 
ent of this model is that the accretion rate decreases with 
decreasing radius from the SMBH, indicating a signifi- 
cant mass loss by a magnetically driven outflow or wind. 
This is necessary to make the model consistent with ob- 
servations, and such wind activities arc also supported 
by recent numerical simulations (YQN03 and references 



therein). 

Though this model can explain the SED of Sgr A* at 
present, there are several lines of evidence for higher ac- 
tivities in the Galactic center in the past, such as much 
higher X-ray luminosity and mass outflows (sec §3 for a 
brief review). Hence, it is interesting to examine whether 
these lines of evidence can consistently be explained by 
changing the accretion rate within the framework of the 
RIAF model. The first aim of this paper is to show that 
it is indeed possible, and I also discuss whether such a 
high accretion rate is reasonable in the environment of 
the Galactic center, as well as the cause of the sudden 
decrease leading to the current low rate. 

The 511 keV electron-positron annihilation line emis- 
sion in the Galaxy has been observed for a long time 
(Knodlseder et al. 2005 and references therein), and the 
latest observation by the INTEGRAL gamma-ray obser- 
vatory revealed that the all sky distribution is dominated 
by the extended bulge component 8° FWHM), with 
a weak evidence for the disk component (Knodlseder et 
al. 2005). Although the disk component can be explained 
by positron emission from radioactive nuclei produced in 
supernovae, the origin of the bulge component is still a 
mystery. A positron production rate of 1.5 x 10"''^ s^^ 
is required to explain the bulge component. 

A number of candidates have been proposed for the ori- 
gin of the bulge component, but the large bulge-to-disk 
ratio excludes many of them related to recent star forma- 
tion activities, leaving type la supernovae (SNe la) and 
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low-mass X-ray binaries (LMXBs) as the primary candi- 
dates (Knodlseder et al. 2005). However, about one order 
of magnitude higher rate of bulge SNe la than the best 
estimate is required to explain the observed annihilation 
rate (Prantzos 2006). The positron production rate from 
LMXBs or microquasars is highly uncertain and might be 
sufhcient (Guessoum et al. 2006), but the bulge-to-disk 
ratio of LMXBs is considerably lower than that of the 511 
keV data. Most of the positrons produced in microquasars 
are probably lost by annihilation before injection into in- 
terstellar medium (ISM) (Guessoum et al. 2006), and 
direct annihilation gamma-ray emission from the sources 
would be inconsistent with the diffuse gamma-ray back- 
ground around MeV band (see §6.3). Annihilation of dark 
matter particles with a mass of ^ MeV is another possi- 
ble solution, but such a particle is not naturally predicted 
by the particle physics theory, in contrast to the well- 
motivated candidates such as the supersymmetric parti- 
cles that are much more massive (m^ ^ 50 GeV, Bertone 
et al. 2004). Furthermore, the upper limit on 511 keV 
line flux from the Sagittarius dwarf galaxy excludes this 
scenario for almost all types of the halo density profile 
(Knodlseder et al. 2005). 

The second aim of this paper is to show that, if the 
past accretion rate onto Sgr A* was in fact higher than 
now with the rate and duration inferred from the obser- 
vational evidence mentioned above, the RIAF model of 
Sgr A* gives a natural explanation for the bulge 511 keV 
line emission. There are models of positron production 
by accretion activity of Sgr A* (Titarchuk & Chardonnet 
2006; Cheng, Chernyshov, & Dogiel 2006), but the model 
proposed here is vastly different from these and it gives a 
more natural explanation for the 511 keV line (§8.1). 

A variety of phenomena on various scales around the 
Galactic center will be discussed in this paper, and these 
are summarized in schematic diagrams in Fig. 1 for the 
reader's convenience. The summary of various comparison 
between observations and the model is given in Table 1. 
Starting from a brief description about the RIAF model of 
Sgr A* (§2), I review and summarize the lines of evidence 
for the past increased activity of Sgr A* , and present RIAF 
interpretations for these, in §3. We discuss the physical 
origin of the past high activity in the Galactic center envi- 
ronment in §4. Then I calculate the positron production 
rate around the event horizon of Sgr A* in §5, and consider 
the ejection from Sgr A* (§6) and propagation in ISM (§7). 
In §8, I compare my model with the other models of the 
511 keV emission, and discuss about future observational 
tests of this model. I apply 8 kpc as the distance to the 
Galactic center (Eisenhauer et al. 2003). 

2. The RIAF Model of Sgr A* and Outflow 

The YQN03 model is a one-dimensional, global RIAF 
model that can reproduce the SED of Sgr A*. Though 
detailed three-dimensional structure is not taken into ac- 
count in this model, models or simulations incorporating 
more complicated structures have not yet succeeded in re- 
producing the Sgr A* SED (e.g., Mineshige, Kusnose, & 



Matsumoto 1995; Ohsuga, Kato, & Mineshige 2005). In 
this paper I discuss the scaling of the radiation and out- 
flow from Sgr A* when the accretion rate is changed, and 
the simple YQN03 model is sufficient and currently the 
best choice for this purpose. 

The accretion rate of the YQN03 model is normalized 
at the outer boundary corresponding to the Bondi radius 
(rs ~ lO^r^ = 0.029 pc), which is inferred from the ambient 
density and temperature of X-ray emitting gas (Baganoff 
et al. 2003), where rg is the Schwartzshild radius of the 
SMBH. The size of the extended (- 1".4 FWHM) X-ray 
emission of Sgr A* (Baganoff et al. 2003) nicely corre- 
sponds to this radius. The accretion rate at this radius is 
Ma^ccirs) ~ ayMs ~ W'^Mq/yt, where ^ 0.1 is the di- 
mensionless viscosity parameter and Mb is the Bondi ac- 
cretion rate. The YQN03 model assumes a radially vary- 
ing accretion rate. Mace oc with s = 0.27, based on the 
adiabatic inflow-outflow solutions (ADIOS, Blandford & 
Begelman 1999) rather than the advection dominated ac- 
cretion flow (ADAF, see Kato et al. 1998; Narayan et al. 
1998 for reviews) that is the simplest RIAF solution with 
a radially constant mass accretion rate. This modification 
from ADAF is necessary since the ADAF model predicts 
too large density and magnetic field strength in the inner- 
most region, in contradiction with the observational upper 
bounds on rotation measure. 

The accretion flow is quasi-spherical and the density 
p is determined by M^^cdf) ~ Anpr'^Vin, where Um is the 
inflow velocity of the accretion flow. In the self-similar 
solution of RIAFs, is close to and proportional to 
the free fall velocity Wff, and hence p oc r~^/^^'^. In 
the YQN03 model, the particle density is rig ~ p/mp ^ 
3.9 X 10''(r/rs)"-'^-^'^ cm~'^, where nip is the proton mass, 
and hence ^ O.lbvs. This also means p oc Mace- 

It is reasonable to assume that the evaporated mass flow 
due to the non-conserving accretion rate is ejected as wind 
or outflow, whose velocity is comparable with the escape 
velocity Uosc at the region where the wind originates. Then 
the mass ejection rate per logarithmic interval of radius is 
=r{dMacc/dr) = sMacc, and the largest kinetic energy 
is produced from the innermost region where Wosc is close 
to the speed of light, c. We denote r ~ for the wind 
production region, and the outflow kinetic energy from 
this region becomes Lkin ~ ^ 1.5 x lO^^r;^"'^^ 

erg/s, where r3 = r^,/ (3rs). 

The jet or outflow may also contribute to the Sgr A* 
SED, especially in the radio bands (see discussion of 
YQN03). The outflow kinetic energy derived here is in 
fact comparable with those assumed in the jet models of 
Sgr A* (Falcke & Biermann 1999; Yuan, Markoff, & Falcke 
2002; Le & Becker 2004), indicating that this estimate is 
a reasonable one. 

The outflow kinetic power is much larger than the X- 
ray luminosity of Sgr A*, which is Lx ~ 10^"^ erg/s in the 
quiescent state and ~ io34-35 gj.gy'g jjj ^j^g flaring state. 
Such a trend of jet-energy dominance is indeed estabfished 
for black holes in stellar X-ray binaries in the low/hard 
state (Gallo et al. 2003, 2005), which is believed to be a 
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state corresponding to RIAF (Esin et al. 1997). In fact, 
the above kinetic and X-ray luminosities are consistent 
with the relation inferred for the stellar black hole bina- 
ries: /kin = /(/jf/0.02)0-5 with / = 0.06-1, where /ki„ and 
Ix are in units of the Eddington luminosity (Gallo et al. 
2005). Though the mass scale of black holes is quite dif- 
ferent, this is expected if (1) the critical accretion rate in 
units of the Eddington rate {M /MEdd) between the stan- 
dard thin disk and RIAF is independent of the black hole 
mass, (2) X-ray luminosity scales roughly as A/f^,^, in the 
RIAF regime, and (3) the mass outflow is proportional to 
Mace ■ The first two are indeed the properties of the AD AF 
(Narayan et al. 1998), and the third is in accordance with 
the radially varying accretion rate assumed in YQN03. 

Some active galactic nuclei (AGNs) show the jet activity 
with much higher speed (the bulk Lorentz factor F ^ 10) 
than the above velocity estimate ('^ Vcsc- only mildly rel- 
ativistic with F ~ 1). However, F ~ 1 is reasonable for Sgr 
A* , since the same trend has been known for the outflows 
from stellar X-ray binaries in the low/hard state, which 
are not strongly beamed and not extremely relativistic 
(F ^ 2), in contrast to those of X-ray transients (Gallo et 
al. 2003; Narayan & McChntock 2005). 

3. Past Higher Activity and RIAF Interpretation 

3.1. Higher X-ray Luminosity 

There are a few independent lines of evidence that 
about 300 years ago Sgr A* was much more luminous 
than now in the X-ray band. Koyama et al. (1996) found 
fluorescent X-ray emission reflected from cold iron atoms 
in the giant molecular cloud Sgr B2 by ASCA observa- 
tion (Fig. 1). Since there is no irradiation source to ex- 
plain the iron line emission, they suggested a possibility 
that ^300 yrs ago Sgr A* was much brighter than now. 
More recent studies by Murakami et al. (2000, 2001a, b) 
found a new X-ray reflection nebula associated with Sgr 
C, and estimated the increased past luminosity of Sgr A* 
as Ljf - 3 X 10^^ erg/s from Sgr B2 and C data. This 
claim was independently confirmed by an INTEGRAL 
observation covering higher energy band of 10-100 kcV 
(Revnivtsev et al. 2004); the ASCA and INTEGRAL 
data of Sgr B2 can nicely be fit by a reflection of radi- 
ation from Sgr A*, whose luminosity is 1.5 x 10"^^ erg/s 
in 2-200 keV and spectrum is a power-law with a photon 
index of /3 = 1.8 ±0.2 {dN/dE (x E-f^). 

"The ionized halo" surrounding Sgr A* with a density 
rie.h ~ 10^-10'^ cm^'^ has been known from radio observa- 
tions (Pedlar et al. 1989; Anantharamaiah ct al. 1999), 
and it extends to a radius of ~ 10 pc (Fig. 1). Macda 
et al. (2002) argued that currently no ionizing source is 
found for the ionized halo, and the past activity of the 
Sgr A* may be responsible for the ionization, requiring 
an X-ray luminosity of Lx ~ 10*° erg/s, a similar flux to 
those inferred from the X-ray reflection nebulae. 

Wc can then estimate the boost factor of the accre- 
tion rate to achieve the X-ray luminosity of Lx ~ 10"^^- 
10*" crg/s in the RIAF model. It should be noted that 
this luminosity is ~ 10~^ in units of the Eddington lumi- 



nosity, and hence it is still well within the ADAF/RIAF 
regime. According to Fig. 5 of YQN04, a boost factor 
fb = Mpast/Mpresont ~ 10^-10* (at & fixed radius r) is re- 
quired to achieve the luminosity of Lx ~ 10^^-10*" erg/s, 
assuming a constant value for s. The SED in YQN04 
with fb ^ 10^ is roughly constant in i/L^, (luminosity per 
logarithmic frequency interval) in the X-ray band, be- 
ing consistent with the X-ray spectrum inferred from the 
X-ray reflection nebula. Then the outflow kinetic lumi- 
nosity in such a higher activity phase should be Lkin ~ 
4.7 X 10*i/3.5r3"°-^' crg/s, where /g.s = /h/lO^ '^. In units 
of the Eddington mass accretion rate, Msdd = lO^Edd/c^ 
(corresponding to a radiative efficiency of 0.1), the accre- 
tion rates at r = rs and r = rg are 4.6 x 10^^/3.5 and 
2.1 X 10"'^ /a. 5, respectively. 

3.2. Mass Outflows on Various Scales 

There is evidence for powerful mass outflow from the 
Galactic center on scales from arcminutes to tens of de- 
grees (Fig. 1). Bland-Hawthorn & Cohen (2003) re- 
ported the mid-infrared emission from dust expanding in 
the Galactic center lobe (GCL) on a few degree scale, and 
estimated the total kinetic energy as ~ 10^^ erg with a 
velocity of ^ 100 km/s and a dynamical time of ^ 10^ yr. 
The size, energy, and time scales are similar to those of 
the expanding molecular ring (EMR) around the Galactic 
center (Kaifu et al. 1972; Scoville 1972). Bland-Hawthorn 
& Cohen further argued that this result is consistent with 
the North Polar Spur (NPS) on an even larger scale (up 
to tens of degrees), which has been interpreted by Sofue 
(2000) to be an outflow from the Galactic center ^, with an 
energy scale of ^ iQ^s-se gj.g ^^^^ dynamical time scale of 
^ 10^ yr. The kinetic luminosity inferred from these vari- 
ous observations is in a nice agreement with the estimate 
for the RIAF model with the boost factor of ~ 10'^^*, 
and hence we can attribute these outflows to the past ac- 
tivity of Sgr A* that was responsible for the higher X-ray 
luminosity. This indicates the duration of ~ 10'' yr for the 
higher activity, which is reasonable compared with various 
estimates for lifetimes of AGNs (Martini 2004). 

The origin of the outflow from the Galactic center may 
instead be starbursts in the nuclear region, as suggested by 
Bland-Hawthorn & Cohen (2003). However, observations 
of gamma-rays from radioactive ^^Al seem to disfavor the 
starburst scenario. The energy of 10^^ erg in 10^ yr is 
equivalent to '--^ 1 supernova per century, almost compara- 
ble with the Galaxy-widc core-collapse supernova rate of 
1.9 ±1.1 per century (Diehl ct al. 2006) estimated from 
the flux of ^^Al gamma-rays, whose spatial distribution 
is clearly associated along with the Galactic disk. Since 
the half-life of ^^Al is 7.2 x 10^ yr, ^^Al cannot travel be- 
yond the GCL region with the outflow velocity of ^ 100 
km/s. Therefore, if the origin of the outflow is starbursts, 
wc expect a strong Al gamma-ray emission concentrated 
within a few degree from the Galactic center, with a flux 

^ There are other interpretations for NPS by closer objects on the 
Galactic disk like a supernova remnant, but see Bland-Hawthorn 
& Cohen (2003) for arguments favoring the Galactic center in- 
terpretation. 
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comparable with the total gamma-ray flux along with the 
disk. However, such a strong concentration at the Galactic 
center is not found (Prantzos & Diehl 1996; Knodlseder 
et al. 1999), indicating that the accretion activity of Sgr 
A* is more plausible as the origin of the mass outflow. 

Muno et al. (2004) studied in detail the Chandra 
data of the diffuse X-ray emission within ~ 20 pc of the 
Galactic center, and they concluded that the hard compo- 
nent plasma {kT ~ 8 keV, Koyama et al. 1989; Yamauchi 
et al. 1990) cannot be explained by unresolved point 
sources. The latest Suzaku observation by Koyama et 
al. (2006) further strengthened the case for the truly dif- 
fuse plasma (but see also Revnivtsev et al. 2006). If it is 
truly diffuse, it cannot be gravitationally bound and re- 
quires a large energy input of ~ 10^" erg/s to keep this 
hot plasma, when the escape time is estimated simply by 
the sound velocity. This is too large to be explained by 
supernova explosions and the origin of this hot plasma 
is still a matter of debate. Later in this paper (§7.1), I 
will discuss in more detail about the consequences of the 
energy injection into ISM by the wind from Sgr A*, and 
show that this hot plasma can be explained as a result of 
shock heating by the wind. 

4. The Origin of the Past Higher Activity 

We have shown that all the lines of evidence for the 
past higher activities of Sgr A* can be explained if the 
mean accretion rate in the past ^ 10^ yrs is higher than 
now by a factor of /& ~ 10'^"^, and Sgr A* had such a high 
rate 300 yrs ago. Then the next questions are: (1) what is 
the source of the accreting matter during the high activity 
phase, and (2) what caused the sudden drop of accretion 
rate by a factor of 10"^"^ on a time scale of just ~ 10^"^ 
yrs. Here I give reasonable explanations for these. 

4-1. The Role of Sgr A East and the Ionized Halo 

Maeda et al. (2002) proposed that, based on their 
Chandra observation of the supernova remnant Sgr A 
East, the past higher activity was induced by accretion 
from the dense supernova shell expanding into the ion- 
ized halo. The location of Sgr A* is in fact inside Sgr A 
East and close to its shell (Fig. 1). In this scenario the 
duration of such high accretion rate is only 10'^ yr, as 
inferred from the shell thickness (~1/10 of the observed 
shell radius r^h = 2.9 pc) and the shell expansion velocity 
{vsh ^ 200 km/s estimated by a simple theoretical model 
of supernova remnants). For comparison, the age estimate 
of Sgr A East is ^ 10'^ yr. A much higher accretion rate 
than now is possible by the Bondi accretion with the shell 
density estimated from shock-compression of the gas of 
the ionized halo. (The supernova ejecta is negligible at 
this stage.) 

However, according to this scenario, we expect a com- 
parable or even higher accretion rate by accretion di- 
rectly from the ionized halo before the passage of the Sgr 
A East shell, since the density enhancement by shock- 
compression is at most by a factor of 4 and the sound ve- 
locity of unshocked gas is likely much lower than the shell 



velocity. Suppose that the SMBH is embedded directly in 
the ionized halo. The Bondi radius in this case could be 
different from the current estimate (0.03 pc) based on the 
X-ray observation (Baganoff ct al. 2003), since the cur- 
rent gas properties around Sgr A* have been altered from 
those of the ionized gas by the passage of the Sgr A East 
shell. The Bondi radius in the halo, rB,h = 2GAjf,/c^^ 
is larger than 1 pc if the temperature kT^ is lower than 
^ 110 eV, where Cs,/t is the sound velocity in the ionized 
halo. A plausible temperature of the ionized halo is 1 
eV (Maeda et al. 2002), and hence the outer boundary of 
the accretion flow would not be determined by the Bondi 
radius, but rather by the gravitational radius rg^av ^ 1 pc, 
within which the gravity of the SMBH is dominant com- 
pared with stars around the SMBH (e.g., Schodcl ct al. 
2003). 

Then we can extrapolate the YQN03 model from tb ~ 
0.029 pc out to Tgrav with the boost factor of ^ 10^"'* 
and rie oc r~^'^^, and the density at rgj-av becomes rie ~ 
1.1 X 10'^ /a. 5 cm^'^. This is consistent with the density 
of the ionized halo, n^^h, and hence the RIAF at the in- 
creased accretion rate is naturally connected to the en- 
vironment surrounding the SMBH. Since this is an ex- 
trapolation of the RIAF solution, the surrounding gas can 
accrete even if it has a significant angular momentum, be- 
cause of the angular momentum loss by viscosity. This is 
in contrast to the simple picture of the spherical Bondi 
accretion. 

In this new scenario, the high accretion rate can last 
for a much longer time scale than ^ 10^ yrs, and accre- 
tion from the dense supernova shell is no longer neces- 
sary. Still, Sgr A East must play an important role to 
explain the sudden drop of the accretion rate ~300 yrs 
ago, by the destruction of accretion flow when the dense 
shell passed through the SMBH. It should also be noted 
that there are the arm-like structures of Sgr A West and 
the circumnuclear disk surrounding it on the scale of ~1 
pc (e.g., Yusef-Zadeh, Meha, & Wardle 2000; see also Fig. 
1), which could be the remnants of the former accretion 
flow. 

4-. 2. Destruction of the Accretion Flow by Sgr A East 

We examine the destruction process more quantitatively 
as follows. Since the ADAF and ADIOS solutions have 
a positive Bernoulli parameter (Narayan et al. 1998; 
Blandford & Begelman 1999), the flow is not gravitation- 
ally bound and change of flow velocity by external force 
would result in a destruction of the flow. Therefore we 
should compare the momentum of the flow and the super- 
nova remnant to estimate the effect of the supernova shell 
passage. The momentum of the accretion flow is: 

Pace ~ 'Y^^P ^flow (1) 

~5.6xl04M^j /3.5gcms-i, (2) 

where we estimated the flow velocity uaow by Vim since 
the rotation velocity is negligible when the adiabatic in- 
dex 7ad ^5/3 in ADAFs (Narayan et al. 1998). The 
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momentum of the supernova remnant given to the flow is: 



where Eg-^ is the shell kinetic energy of the supernova 
remnant. Therefore the accretion flow could have been de- 
stroyed at r ^ 0.1 pc from the SMBH. Destruction should 
have occurred with a time scale of the shell crossing (~ 10'^ 
yrs), and the accretion time scale at this radius is also 
r/vin ~ 1.3 X 10"^ yrs. These time scales are consistent 
with the required time scale of the accretion rate drop, 
^ 300 yrs. The current low rate may be determined by ei- 
ther the diffuse gas in the supernova remnant, winds from 
nearby stars (e.g., Cuadra et al. 2006), or residual of the 
former accretion flow. 

4-3. Comparison with Nearby Galaxies 

We may ask how the suggested high activity of Sgr 
A* in the past compares with those found in nearby nor- 
mal galaxies, because it would be statistically unlikely if 
our Galaxy had a much higher activity compared with 
nearby normal galaxies. As reviewed by Ho (2004), nu- 
clear activity is quite commonly found in nearby galax- 
ies. The YQN model with /{, ^ 10'^ predicts a nuclear 
B-band luminosity of ~ 10"^^ erg/s, and the number 
density of such galaxies expected from the nuclear lu- 
minosity function is ~ 10~^(/)./0.75)'^ Mpc""^, which is 
similar to that of typical galaxies like our own, where 
h ~ i7o/(100 km/s/Mpc) is the Hubble constant. Nuclear 
radio cores with flux of 10^^-10^^ W Hz~^ at 5 GHz are 
also commonly found, which is again a similar radio lu- 
minosity to that predicted by the YQN04 model with 
fb ~ 10'^"'^. Therefore, the increased activity of Sgr A* is 
not particularly rare compared with nearby normal galax- 
ies, indicating that the characteristic time scale of the in- 
creased activity can be 10^ yr or even longer, possibly as 
long as the cosmological time scale. The accretion rate at 
the event horizon, i.e., the mass growth rate of the SMBH 
is 1.4 X lO~'*/3.5M0yr~^, and the mass growth in 10 Gyr 
is 1.4 X 1O^/3.5M0, which is still less than M,. 

AGN activity is generally sporadic and showing strong 
variability, and hence it is naturally expected that the ac- 
cretion rate was changing significantly always in the past 
^ 10'' yrs. In fact, the characteristic structures such as 
GCL, EMR. and NPS indicate such variability or modu- 
lation of the accretion and wind activity. However, con- 
clusions of this paper are not seriously affected if the mean 
or characteristic accretion rate is given by fi, ~ 10'^^"'. 

4-.4- VVe Living in a Very Special Time? 

One may feel uneasy if we are witnessing a large drop 
of the accretion rate just during ^300 yrs compared with 
the normal accretion time scale of ~ 10^ yrs. indicating an 
extremely low chance probability of observing Sgr A* in 
the present phase: 300/10'' = 3 x 10~^. Here I argue that 
the actual chance probability is likely to be much larger 



than this simple estimate. 

About 100 young massive stars arc known within the 
central parsec of Sgr A* (Paumard et al. 2006). Some 
of them are already in the post-main-sequence phase, and 
the age of this stellar population is estimated as ^ 6 Myr. 
This indicates a supernova rate of ~ 1.7 x 10~^yr~^, and 
the chance probability of observing Sgr A* within 300 yrs 
after the shell passage becomes 0.5%. This is small, but 
not extraordinarily small for realization. The mean ac- 
cretion rate can be kept at the level of /f, ^ lO'^"'* if the 
accretion rate comes back to the mean level within the 
typical time interval of supernovae, ~ 6 x 10^ yrs, after 
the shell passage. 

It is an observational fact that we are living in a some- 
what special time, if the association of Sgr A* with the 
Sgr A East shell is not just a superposition in the sky but 
physical. The age of ~ 10^ yrs of Sgr A East indicates 
that the supernova rate within the central parsec of Sgr 
A* should be less than 10~^ yr~^- Therefore, the chance 
probability must be less than 300/10'* = 3%, which is just 
6 times larger than the above estimate. Therefore I con- 
sider that the argument of the chance probability does not 
seriously weaken the case for the overall picture proposed 
by this work. 

5. Positron Production around Sgr A* 

5.1. Physical Quantities around the Event Horizon 

Now I consider the pair-production at the region where 
the wind originates, r ^ r*. During the increased phase, 
the particle accretion rate is iVacc = -^acc/^p ~ 7.1 x 
lO'^'^/a.sfg'^^ s~* and the particle outflow rate per loga- 
rithmic radius is = siVacc- The particle density in the 
accretion flow is rig ~ 3.2 x 10-^°/3.5r;^*'^^ cm~^. The ac- 
cretion time spent around this radius is iacc T*/'Uin ^ 
1.0 X lO^rg^^ s. 

The electron temperature of the YQN03 model is Tg ~ 
8 X lO^^rg"* K, and hence relative motion of electrons 
is sufficiently relativistic at the pair-production region 
if ~ 3rs, with the mean electron Lorentz factor ~ 
3.151A:Te/(TOeC^) ^ 40r^* in the rest frame of the accre- 
tion flow. The temperature does not heavily depend on 
the enhancement factor /j,^ 10^'^, if the transfer efficiency 
of viscous heating energy from ions to electrons (the pa- 
rameter 6 in YQN03) does not change with the accre- 
tion rate, as assumed by YQN04. It may increase with 
increasing accretion rate because of higher density and 
hence more efficient interactions, but the value of S as- 
sumed in YQN03 is already of order unity (=0.55), not 
leaving much room for the increase of S. 

5.2. Pair Equilibrium Criterion 

The positron density produced in the accretion flow de- 
pends on whether the pair production process achieves 
the equilibrium with the pair annihilation. This 
can be evaluated by comparing the pair production rate 
density and pair annihilation rate density ?i±,ann = 
nen_|_cr±_annC, whcrc rt-|_ is the produced positron density. 
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The pair annihilation cross section is 

9 

7rr 

(T±,ann = — [ln27±-l] (5) 

7± 

in the ultra-relativistic limit (Svcnsson 1982), where is 
the classical electron radius and 7± is the Lorentz factor 
in the rest frame of one particle. Hence 7± can be related 
as 7± ~ 7+7e, where 7+ is the positron Lorentz factor in 
the flow frame. Initially the positron density is small, and 
it increases with time as ^ n^t until t ~ t±,ann when 
the equilibrium is achieved (n+ = J^i.ann), where the pair 
annihilation time scale is t±,ann = {neca±^s.nn)~^ ■ Since 
^ifc.ann docs not dcpcud on the equilibrium condition 
is the same for any pair production processes and it is de- 
termined by comparing t±,ann to the accretion time scale, 
as: 

= 4.2 X 10-2^- V3.5rl■^^ (g) 

lb, ann 

where we estimated the logarithmic part of the cross sec- 
tion by 7+ ~ 7e ^ 40. 

For positrons produced by electron-electron scatter- 
ing (e^e^ e^e^e~e+), we expect 7+ ~ 7e, while for 
positrons produced by two photon annihilation (77 — > 
e~e^), a variety of 7+ is possible (see the following subsec- 
tions). The positron energy may also significantly change 
by interaction with the accreting material within the ac- 
cretion time scale. However, for any value of 7+, we find 
that <acc ^ i±,ann, and hence the equilibrium will not be 
achieved, meaning that we can estimate the pair density 
by n+ 7i+tacc- 

5.3. Spectral Energy Distribution of Sgr A* 

For pair production processes including photons, we 
must assume the form of the SED and the emission re- 
gion of Sgr A* during the increased activity phase. Here 
I assume that the radiation is mainly coming from the re- 
gion around the event horizon. It should be noted that 
the X-ray emission in the quiescent phase of the YQN03 
model for present-day Sgr A* is dominated by thermal 
bremsstrahlung at large radii far from the SMBH, which 
is in agreement with the extended X-ray emission (~ 1"). 
However, the synchrotron-self-Compton (SSC) component 
becomes dominant when the accretion rate is increased as 
fb ^ 10'^ (YQN04), and it is produced in the innermost 
region. Therefore the above assumption is reasonable for 
the higher activity phase. 

For calculations below, I assume a constant SED in i^L,y. 
The predicted SED of the YQN04 model when > 10^ 
is approximately flat in i/Li, from keV to MeV band. A 
flat i/Lij SED is also supported in 2-200 keV band by the 
observed spectrum of the X-ray reflection nebula (§3.1). 
Therefore, this assumption is well supported both by ob- 
servation and theory in the photon energy band of keV- 
MeV. 

The SED beyond MeV is more uncertain, and here I 
simply examine the constraints on the current SED of Sgr 
A*. The TeV gamma-ray emission detected by H.E.S.S. 
from the inner 10' of Sgr A* (Aharonian et al. 2004) is 



possibly coming from the region close to the event hori- 
zon, and its flux is similar to that extrapolated from X-ray 
bands with a constant vL^, spectrum. Though the GeV 
flux of 3EG J 1746- 2851, which is the closest to Sgr A* 
among the EGRET sources, is considerably higher than 
the X-ray flux of Sgr A*, the poor angular resolution in 
this band does not allow to establish a firm connection 
between the GeV flux and Sgr A* (Mayer-Hasselwander 
et al. 1998; Aharonian & Ncronov 2005)^. Therefore, it is 
not unreasonable to assume a flat uL^, SED at photon en- 
ergy of ^ MeV, though uncertainty is large. On the other 
hand, the photon production region will become optically 
thick to pair production for very high energy photons 
(^ GeV), and the constant vLi^ assumption will be no 
longer valid at such high energy band (see §5.4.3). 

5.4. Expected Pair Amount in the Wind 

Now I estimate the expected amount of pairs pro- 
duced by three processes of pair production, i.e., electron- 
electron scattering {e~e~ —^ e~e~e~e~^), photon-electron 
collisions (76^ e^e^e^), and two photon annihilation 



(77- 



The rates of corresponding proton processes 



(e.g., p e ^p e e e+) are about one order of magnitude 
smaller than these (Zdziarski 1982, 1985). 

5.4.1. Electron- Electron Scattering 

I used the formula given in Svensson (1982) for the 
cross section in the ultra-relativistic limit, which is aee = 
1.7 X 10~^* cm^ for 7e = 40 and depends on 7e only log- 
arithmically. Hence I ignore the dependence on j^. Then 
the density ratio of the produced positrons to electrons is 
given as: 



n+ 

Up, 



n+ta 



Ca-eeTlet 



He 



1.6 X 10-V3.5r^ 



0.27 



(7) 
(8) 



and hence the total positron production rate as an outflow 
from Sgr A* is: 



= 3.2 X lO'^Vs 



41 f2 ^0.54 



5'3 



(9) 

(10) 

This is a rate per Inr*, and integrating from = 3rs to 
40 X 3rs, beyond which electrons become non-relativistic 
and hence the above formulations arc no longer valid, 
the rate is increased by a factor of 11.7 leading to iV+ ^ 
3.7 X 10^^/1 5 s~^. This estimate is, within the model un- 
certainties, in nice agreement with the rate required for 
the bulge 511 keV emission, 1.5 X lO''^ (Knodlseder et 
al. 2005). 

5.4.2. Photon- Electron Collisions 

The cross section for e~j — > e~e~e^ depends on the 
photon frequency Vci measured in the electron's rest 

^ A recent analysis of the EGRET data by Hooper & Dingus 
(2004) excluded Sgr A* as the origin of 3EG J1746-2851 beyond 
the 99.9% confidence level. 
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frame, which is given by (Svensson 1982)^: 
in the non-rclativistic limit (xcr — 4 <C 4), and 



(18) 



28 

— ln{2xc 



218 
"27 



(11) 



(12) 



in the ultra-relativistic limit (xor — 4 3> 4) , where x^r = 
hvcr / (TUeC^) , a is the fine structure constant, and the re- 
action threshold is Xd-.th = 4. 

Treating electrons as a single energy population with 
7e ~ 40, and estimating the photon number density per 
unit photon frequency [v) in the laboratory frame as 

^ / {'^■nrlchv), the positron production rate by this 
process can be written as: 



ycdv (13) 

■ 22/ 2 '^^Jer , (14) 

where we have used ^^v. With the assumption of 

a flat vL^ SED. the integration over Xa- is mostly con- 
tributed from photons with Xa ~ 20, i.e., hi/ ~ 0.3 MeV 
for 7e ~ 40. This target photon energy is within the range 
where the Sgr A* luminosity during the increased activity 
can reliably be assumed, and hence the uncertainty about 
the luminosity and SED is small. 

Then the produced positron density is given by: 



Up 



n+t. 



+ ^acc 



4.2 X 10"%9.5r-^-^ 



(15) 



where L39.5 = iyL^/{10 erg/s), and the total positron 
production rate in the outflow is: 



A^4 



8.1 X 104%9.5/3.5r3-l-23 g-1 



(16) 
(17) 



Integration over Inr* at > Sr^ would slightly decrease 
the above number by a factor of 1/1.23. Though this 
number is smaller than the rate required to explain the 511 
keV emission by about one order of magnitude, it may also 
be important if the model uncertainties are considered. 
Note that oc M^^^ at photon energies higher than the 
X-ray band in RIAFs, and hence cx 



5.4-3. Photon-Photon Annihilation 

The pair-production by two photon annihilation most 



efflciently occurs with a cross section of cr-y 



1.7 X 



10^^^ cm~^, when the photon energy at the center-of- 
mass is about the electron mass, as {hvihvhY/'^ ^ 2meC^, 
where vi and Vh are the frequencies of two photons at the 
laboratory frame [vi < Vh) (e.g., Salamon & Stecker 1998). 
For photons meeting this condition, the pair-production 
rate density is given by 

^ The numerical factor Sy^ of the non-rclativistic formula in 
Svensson (1982, eq. 31) should be corrected to ttVS as noted in 
Svensson (1984). 



= 3.4 X 10"%9.5r3"^ 



which is constant against i>h (and ui) by the assumption of 
the constant i^L^ SED. It should be noted that this is valid 
only when the region is optically thin for high frequency 
photons, i.e., r^^ ^ 1, where 

(19) 

, . (20) 
IMeV/ ^ ' 

The luminosity and pair production would then be sup- 
pressed for very high energy photons of GeV.* Therefore 
we expect that the pair production rate will mostly be 
contributed by photons in keV-GeV bands. 
Now the density of positrons produced is: 

Tie Ue 

= 8.4xlO-^Ll,,J^lrj'-'\ (22) 

and hence the total positron production rate from the Sgr 
A* is: 



Hp 



l.eXlO^lL^g.grg-l S"! 



(23) 
(24) 

Note that this estimate is per unit logarithmic interval 
of i/h (or i^i). Since the Sgr A* luminosity during the 
increased activity can reliably be modeled up to ^ MeV, 
the uncertainty is small for photons of hvi ~ hiy^ ~ mgC^ . 
Though it suffers larger uncertainty about the SED in 
MeV-GeV bands, integration over hiyfi ~ 1 MeV-1 GeV 
would increase the total rate by a factor of 6.9. The rate 
has a large dependence on the accretion rate as cx 
(X f^. Considering the model uncertainties, the two 
photon annihilation may also substantially contribute to 
the observed bulge 511 keV emission. 

5.5. Comparison with Previous Studies 

Positron production in RIAFs has been investigated 
in several previous studies (Kusunose & Mineshige 1996; 
Bjornsson et al. 1996; Esin et al. 1997), and the re- 
sults presented here are consistent with these in a sense 
that the effect of pair production in RIAFs is small, i.e., 
n+ ^ Up. However, the values of (n+/ne) found in this 
work are quantitatively higher than the results of Esin et 
al. (1997) for stellar X-ray binaries, who found that the 
pair annihilation is in equilibrium with the pair produc- 
tion by the electron-electron collisions, with the density 
ratio of n^/ue'S 10^^. The difference can be understood 
as follows. 

As discussed in §5.2, the criterion of the equilibrium 
between pair production and annihilation is given by 
tacc/t±^e,nn- According to thc RIAF theory, I find the 
dependence of this quantity on the relevant parameters 
as: 



Because of the increased luminosity compared with that of Sgr 
A* at present, this energy scale is much smaller than the esti- 
mate by Aharonian & Neronov (2005) for the present-day Sgr 
A* (~ 10 TeV). 



8 



TOTANI 



M 



M 



OC 0-±.ann -jy °^ ^i.ann , ^ 



(25) 



Edd 



Though the Eddington accretion ratio AI/M-Edd consid- 
ered in this work is within the parameter range investi- 
gated by Esin et al. (1997), c±,ann IS Considerably dif- 
ferent by the difference of the electron temperature. In 
the two-temperature ADAF model considered by Esin et 
al., the electron temperature is about the virial value as 
Te - 10^-^° K, while it is much higher in the RIAF model 
of Sgr A*, as Te ^ 10^^ K, because of the high value of S 

(§5.1). Since Cr±,ann OC 7^:^ 7~2, the value of tacc/i±.ann 

for ADAF is much higher than the model presented here, 
leading to the pair equilibrium (tacc > i±,ann) as found by 
Esin et al. 

In the pair equilibrium, the positron density is given by 
^ n_i_/(necr±,annc). In the case of the pair production 
by electron-electron collisions (ri+ ~ n^CTeec), we find: 



n+ 



0'±.a 



(26) 



In the model presented here, this value is n^/ue ^ 0.15 
for 7e ^^40, and hence n^/ne> 10~^ is possible even if the 
pair equilibrium is not reached. On the other hand, (T±^ann 
is much larger and a^e is smaller for the two-temperature 
ADAF having a lower Te, leading to the low saturation 
value of (n+/ne)oq ^ 10~^ as found by Esin et al. 

Therefore, the higher electron temperature of the RIAF 
model of Sgr A* than the ADAF model, which is inferred 
from the fitting to observations, is essential to achieve 
the high positron production rate to explain the 511 keV 
emission. 

6. Positron Ejection from Sgr A* 

6.1. Annihilation at the Production Site 

Some of the positrons are lost by annihilation around 
the production site, producing direct annihilation gamma- 
ray emission from Sgr A*. Its spectrum depends on 
the positron spectrum and gravitational redshift, and 
it is thermal with a temperature Tg ~ 10 MeV for the 
electron-electron scattering. However, the small value 

of (tacc/t±,ann) ~ "-±,ann/'^+ ^ 4.2 X 10~^7+^ indicates 

that most of the produced positrons are conveyed into 
the SMBH or ejected by the wind. The direct annihi- 
lation gamma-ray luminosity from the pair production 
site, L±^ann ~ (4t^/3) n±.ann, IS tlicu mucli smaller than 
the positron production and ejection rate by the wind, 
{n+/ne)Nw ^ 47rsrij7i+. 

6.2. Annihilation in the Wind 

Here we check that the positrons trapped in the outflow 
are not lost by annihilation before escaping the SMBH 
gravity. We assume a constant wind speed as ~ Vescir*) 
and the density of the outflow is ~ Nw/{4:7^r^Vw) for a 
steady wind. Then the fraction of positrons that are lost 
by annihilation during wind propagation is: 

r dr 

f/ann = / ?^toO■±,annf ± , (27) 



where v± is the relative velocity of electrons and positrons. 
Initially both electrons and positrons are relativistic, and 
hence v± '--j c and (T±,ann oc 7_|r^ ~ (7e7+)~^ c< r^/'^, since 

7e and 7+ scale as oc n]J,^ cx r~^/^ by adiabatic expan- 
sion. Either electrons or positrons become non-relativistic 
at Tnr = r* min(7e*, 7+*)'^/^, and the scaling changes as 
Ci.ann OC T^^^ at r > Tnr, whcrc 7e* and 7+* are the Lorentz 
factor at the wind creation site (r r»). Then, cr±^annW± 
becomes constant after both electrons and positrons be- 
come non-relativistic, since (T±,ann = t^t^I (v±l c) in the 
non-relativistic limit (Svensson 1982). Therefore, the 
main contribution to the integration comes from r ^ r,ir, 
and we obtain: 

^ann ~ 6ri,u,(r„r)cr±,ann(?'nr)7 TTT^'ni- (^8) 



' 6n^„(r,)cr±^ann(r*)- 



1/3 



(29) 



1.5 X 10-3/3,5-^_ 



+ *'3 



[min(7eH.,7+*)]^^^ • (30) 



This is sufficiently small for any value of 7+*, and hence 
we expect that almost all the positrons produced around 
the SMBH will escape from the SMBH gravity field once 
they are trapped in the outflowing material. 



6.3. On the Constraint from 
Background Radiation 



the Gamma-Ray 



The above results indicate that the annihilation lumi- 
nosity directly from Sgr A* is expected to be much smaller 
than the bulge 511 keV line emission, in the phase of the 
past higher activity. The direct annihilation luminosity at 
present is even much smaller by the boost factor /t,, far 
below the detection limit of gamma-ray telescopes in the 
foreseeable future. 

This is important concerning the constraint from 
the diffuse gamma-ray background, whose flux is 
E {dF/dE) ~ lO"'^ photons cm'^s"^ at -1 MeV within 5° 
from the Galactic center (Beacom & Yiiksel 2005). If any 
sources of the positrons directly emit annihilation gamma- 
rays before injection of positrons into ISM, their spectrum 
would be broad around ~ MeV reflecting the positron 
spectrum and gravitational redshift. The flux from this 
process should not violate the observed MeV background, 
leading to an upper bound on the annihilation rate as 
^ 3.9 X 10^^ s~^, assuming no positronium formation for 
annihilation within the source. Comparing to the positron 
annihilation rate of ~ 1.5 x lO^'^s"^ in ISM, the annihila- 
tion rate within the sources must be ^ 2.6% of the annihi- 
lation rate in ISM, if the sources are concentrated within 
5° of the Galactic center. If the sources are distributed 
like the observed 511 keV photons, the constraint becomes 
< 10%, taking into account that --24% of all the 511 keV 
photons come from the region within 5°. 

The model presented here well satisfies this constraint, 
because of the current low accretion rate. However, it 
puts a stringent constraint on another explanation for the 
511 keV emission by accreting black holes, i.e., LMXBs or 
microquasars. Theoretical models of pair production in 
these objects predict that most (~90%) of the produced 



No. ] 

pairs annihilate near the production site before injection 
into ISM (Guessoum et al. 2006), which is in serious con- 
flict with the above constraint. 

7. Positron Propagation in Interstellar Medium 

7.1. Dynamics and Energetics of the Wind Injected into 
ISM 

When positrons escape from the gravitational poten- 
tial well of the SMBH, the kinetic energy of the outflow 
is expected to be dominant compared with the thermal 
energy because of the adiabatic cooling, and hence the 
positron energy is determined by the bulk Lorentz fac- 
tor r ~ 1 of the wind, i.e., ~ 1 MeV, as argued in §2. 
This is important, since too relativistic outflow would be 
inconsistent with the observational upper bound on the 
injection energy, < 3 MeV (Beacom & Yiiksel 2005; Sizun 
et al. 2006). This bound has been derived by a similar 
argument in the previous section (§6.3), i.e., by requiring 
that gamma-rays from in-flight annihilation of positrons 
before thermalization in ISM do not make an excess of the 
diffuse gamma-ray background radiation. 

The wind will sweep up ISM and heat it up by shocks. 
The ram pressure of the wind will be balanced with ISM 
at Tp, which is determined as N.^nipV^, / (Attt^) ^ Pjsmi as- 
suming a quasi-isotropic wind. Assuming Pism = -S^/(87r) 
by interstellar magnetic field oi B ^ 10/iG in the bulge 
region (LaRosa et al. 2005; Jean et al. 2006), we find 
Tp ~ 3.4 X IQ^fl'^r^^ '^'^B^^ pc, where = B/(10/iG). 
This radius is comparable with the size of GCL or EMR, 
and also with the FWHM of the spatial extent of the 511 
keV emission. Beyond r ~ Tp, the shock-heated gas will 
expand by thermal pressure. Therefore, even if the wind 
originally has some anisotropy, it will not directly appear 
on a scale larger than ^ Tp. 

The observed expansion velocity of ~100 km/s at the 
GCL/EMR region indicates an escape time of tesc ^ 10® 
yrs from this region. Then, the wind kinetic energy stored 
in the GCL region is 

-BqcL ^ -^kin^csc (31) 

-1.5xl0-/3..3-»-(j^)»g. (32) 

This is interestingly similar to the energy of the hot gas 
[kT ^ 8 keV) in the Galactic center observed by X-rays 
(Koyama et al. 1989; Yamauchi et al. 1990; Muno et 
al. 2004), E^ot 2.6 x 105'*(rp/300 pc)^/^ gj.g^ ^^^^^^ j 
obtained this value from surface energy density estimated 
by Muno et al. (2004) assuming the size and depth to 
be ^ Tp. In fact, the observed size of the hot gas 1.8° 
FWHM, Koyama et al. 1989; Yamauchi et al. 1990) is 
in good agreement with Tp. Therefore, the large amount 
of energy stored in the hard X-ray emitting gas can be 
explained by the wind activity. The typical cooling time of 
the hot gas is ~ 10® yr (Muno et al. 2004), which is much 
longer than iosc, and hence the expansion is adiabatic as 
argued by Sofue (2000) and Bland-Hawthorn & Cohen 
(2003). The expansion velocity (~ 100 km/s) is much 
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lower than the sound velocity of the hot gas 10'^ km/s), 
but it is possible if the associated cold material works 
as a ballast, as inferred from the infrared emission from 
expanding dust (§3.2). 

7.2. Positron Propagation Distance 

At r^Tp, positrons are expected to interact with ISM 
and produce the 511 keV emission. In the hot phase of 
ISM, which is the dominant component in the volume fill- 
ing factor in the bulge region (Jean et al. 2006), positrons 
are thermalized in a time scale of 3 x 10® yrs and then 
annihilate in a time scale of ~ 10^ yrs (Guessoum et al. 
2005; Jean et al. 2006). Therefore, the time scale of the 
past higher activity inferred from the outflow evidence is 
sufficient to explain the 511 keV line emission. Positrons 
can maximally reach ~ 1 kpc by the large scale outflow 
(~ 100 km/s) within this time scale, being consistent with 
the observed maximum extent of the 511 keV emission 
(~ 20°). Most positrons must annihilate before traveling 
this distance, since the spectral analysis of the 511 keV 
line indicates that almost all positrons are annihilating in 
warm neutral or warm ionized phase of ISM, where the an- 
nihilation time scale is much shorter. Exact spatial profile 
of the 511 keV line emission would be determined by the 
probability of positrons entering into the warm phase of 
ISM at r > Tp. 

Diffusion in random magnetic field should also have sig- 
nificant effect. Jean et al. (2006) estimated the propaga- 
tion length by quasilinear diffusion for MeV positrons in 
the bulge as ^ 260 pc for a time scale of ^ 3 x 10® yrs, using 
a diffusion coefficient D ^ 9.8 x 10""* kpc^Myr"* that was 
derived from _B = 10 /iC and the Kolmogorov turbulent 
spectrum. Though the Kolmogorov spectrum is not valid 
for cosmic ray propagation in the Galactic disk (Maurin 
et al. 2001), this is not significant because a similar dif- 
fusion coefficient of 4.1 x 10"''' kpc^Myr"-'^ is obtained by 
using the parameters derived by Maurin et al. (2001) to 
explain the cosmic ray data around the Solar system. 

A more detailed, quantitative prediction about the 511 
keV line morphology is beyond the scope of this paper, 
but these considerations indicate that the model presented 
here is broadly consistent with the observed morphology 
and spatial extent of the 511 keV line emission. 

8. Discussion 

8.1. Comparison with Other Models of the Positron 
Production from Sgr A* 

Titarchuk & Chardonnet (2006) proposed a scenario in 
which positrons are produced by annihilation of hard X- 
ray and ^ 10 MeV photons around Sgr A*. The hard 
X-ray photons are emitted from the accretion activity of 
the SMBH, while the ~ 10 MeV photons are produced 
by accretion onto hypothetical small mass black holes 
(SmMBHs) with a mass ~ 10^^ g (- IO-^^Mq), which 
are assumed to be distributed within r ~ 10^~"^rs of the 
SMBH. If they are accreting with the Eddington accretion 
rate and their density is comparable with the dark mat- 
ter, the SmMBHs can supply enough 10 MeV photons for 
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the required pair production. Such SmMBHs have not yet 
been excluded as a candidate of the dark matter, but there 
is few theoretical support in contrast to the well-motivated 
dark matter candidates such as neutralinos (e.g., Bertone 
et al. 2004). Furthermore, it is extremely difficult to 
supply the accreting material onto such SmMBHs at the 
Eddington rate by Bondi accretion; the Bondi accretion 
rate for a 10^^ g SmMBH with typical parameters are: 

a:/.^7.9x10-(£^)(^) ''%/s, (33) 

which should be compared with the Eddington accretion 
rate of MEdd ~ 72 g/s. The difference is almost 24 or- 
ders of magnitude, and it seems quite unlikely that such 
SmMBHs have the Eddington accretion rate at any real- 
istic astrophysical circumstances. 

Cheng, Chernyshov, & Dogiel (2006) considered cosmic- 
ray production by jet or outflow from the SMBH, which 
is ejected when stars are captured by the SMBH. In their 
scenario, the cosmic rays produce pions by collisions with 
ISM, and then positrons are produced by pion decays. 
The most important difference of this model from that 
proposed in this paper is the high injection energy into 
ISM (^ 30 MeV) of pion-decay positrons. Such a high in- 
jection energy is inconsistent with the observational upper 
bound on the injection energy, ^ 3 MeV, as mentioned in 
§7.1. 

Another problem about positron production from pion 
decays is the observed large bulge-to-disk ratio of the 511 
keV line emission. We know that the Galactic gamma-ray 
background in the GeV band is mainly composed of pion- 
decay gamma-rays produced by cosmic-ray interactions in 
ISM (Strong et al. 2000, 2004), i.e., the same process with 
the model of Cheng et al. The GeV background is clearly 
associated along with the Galactic disk, and we have a 
difficulty to explain why we do not see the strong disk 
component of the 511 keV emission if it is produced by 
the cosmic-ray interactions. 

8.2. Predictions and Possible Tests by Future 
Observations 

Although a quantitative prediction of the morphology 
of the 511 keV emission is beyond the scope of this pa- 
per, it could be more spherically asymmetric compared 
with other explanations such as SNe la or MeV-mass dark 
matter. Asymmetry is expected by the wind anisotropy 
in the region of r ^ fp ~ a few degree, or by matter dis- 
tribution in the Galaxy at r ^ rp. Asymmetry has not 
yet been detected in the observed 511 keV morphology 
(Knodlseder et al. 2005), but it does not reject the model 
presented here because of the limited angular resolution of 
the SPI spectrometer of INTEGRAL (- 3° FWHM) and 
the large uncertainty of theoretical expectation for the 
anisotropy. Future observational studies on the 511 keV 
morphology with better angular resolution might, how- 
ever, detect larger asymmetry than that expected for the 
other explanations. 

Another prediction of the proposed scenario is that it is 



extremely difficult to detect 511 keV emission in regions 
other than the Galactic bulge from the source population 
that is responsible for the bulge component. SMBHs are 
generally found in galaxies having bulges, and the near- 
est SMBH is probably M31* in the Andromeda galaxy. 
The X-ray luminosity of M31* is ~ 10^^ erg/s (Garcia et 
al. 2005), which is about 10^ times larger than the quies- 
cent phase of Sgr A*, but 10^ times smaller than the past 
higher activity phase. If this X-ray luminosity reflects the 
typical activity averaged over a time scale of ^ 10^ yr, we 
expect that the 511 keV luminosity of the M31 bulge is 
much fainter than that of the Galaxy. The 511 keV emis- 
sion has been detected at ~ 50ct level (Knodlseder et al. 
2005), and considering the distance to M31 (770 kpc), a 
large improvement of the sensitivity is required. 

On the other hand, there is a better chance to detect 511 
keV emission from regions other than the Galactic bulge, 
for some other models of the bulge 511 keV emission (e.g., 
Knodlseder et al. 2005). For example, if the origin is SNe 
la, improved instruments in the near future will detect 
511 keV line emission from nearby supernova remnants. 
In fact, an interesting limit on the positron escape frac- 
tion from SN 1006 has already been obtained (Kalemci 
et al. 2006). If the origin is the MeV-mass dark matter 
annihilation, we expect 511 keV line emission from nearby 
dwarf galaxies by a modest improvement of the sensitivity. 
If future negative results ruled out the other explanations 
of the 511 keV emission, it would strengthen the case for 
the scenario presented here. 

It should also be noted that the standard prediction 
for positron production from SNe la is only one order of 
magnitude short of that required to explain the bulge 511 
keV emission (Prantzos 2006). Therefore, just a detection 
of 511 keV emission from a nearby SN la does not confirm 
SNe la as the origin of the bulge component, but a close 
examination of the positron production efficiency would 
be required. 

9. Conclusion 

In this paper, I have shown that the several independent 
lines of evidence for the past higher activity of the Galactic 
center (i.e., higher X-ray luminosity and large scale out- 
flows) can quantitatively be explained by the RIAF model 
of Sgr A*, in which energetic outflow plays an essential 
role. A single increased accretion rate from the current 
value explains both the past high X-ray luminosity and 
kinetic luminosity of outflow inferred from observations. 
The required boost factor of the accretion rate is about 
-|^q3-4 ^qj. ^ tijjjg scale of ~ 10^ yrs in the past. I have 
shown that this accretion rate and its duration are natu- 
rally expected in the environment of the Galactic center. 
The outflow is energetic enough to supply the heat to keep 
the hot ^ 8 keV plasma observed in the Galactic center 
region, for which there have been no good explanation. 
The current low accretion rate is likely caused by a sud- 
den destruction of the accretion flow when the shell of the 
supernova remnant Sgr A East passed through the SMBH 
about 300 yrs ago. The chance probability of observing 
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Sgr A* in such a low activity phase is estimated to be 
~ 0.5%. 

I then estimated the production rate of positrons during 
the high activity phase, which are created in the region 
around the event horizon and then ejected by the outflow. 
The rate is found to be comparable with that required 
to explain the 511 keV line emission toward the Galactic 
bulge. I considered three processes of pair-production 
via electron-electron scattering, photon-electron collision, 
and photon-photon annihilation, and interestingly all the 
three processes give a similar positron production rate. 

Therefore the model presented here gives a new and 
natural explanation for the 511 kcV line emission toward 
the Galactic bulge. The favorable aspects of this model 
are: (1) the correct positron production rate, (2) the large 
bulge-to-disk ratio and correct spatial extent in the bulge, 
(3) (i) negligible annihilation near the positron production 
site before injection into ISM, and (ii) the injection energy 
of ~ MeV, both of which are important to meet the con- 
straint from the MeV gamma-ray background, and (4) no 
exotic assumptions or parameters. To my knowledge, the 
other explanations for the 511 keV emission proposed so 
far do not satisfy all of these. 

Anisotropy of the morphology of the 511 keV emission 
larger than that expected for other models would be a 
signature for the model proposed here, which might be 
revealed by future observations with better angular reso- 
lutions. Detection of 511 keV lines from centers of nearby 
galaxies by the same mechanism will be difficult even in 
the foreseeable future. In contrast, some other scenarios 
predict detectable 511 keV lines in directions other than 
the Galactic center by a modest improvement of the sen- 
sitivity. 

I would like to thank an anonymous referee for useful 
comments. This work was supported by the Grant-in- 
Aid for the 21st Century COE "Center for Diversity and 
Universality in Physics" from the Ministry of Education, 
Culture, Sports, Science and Technology (MEXT) of 
Japan. 
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Table 1. Summary of observations versus the model 






Observation Scale (pc) Requirement* 


Prediction^ 


Section 


Evidence for higher past X-ray luminosity 


X-ray reflection nebulae 100 Lx ~ 3 x lO'^'^' erg/s 
Ionized halo 10 Lx ^ 10^° erg/s 


5 X 1039/3.5 
5 X 1039/3.5 


3.1 
3.1 


Evidence for past outflow activity 


Galactic center lobe (GCL) 200 Lkin ~ 3 x 10^^ erg/s 
Expanding molecular ring (EMR) 300 Lkin ~ 1 x lO^^ erg/s 
North polar spur (NPS) 4 x 10^ ikin 1 x lO^i erg/s 
8 keV diffuse gas 300 Ei^at ^ 2.6 x 10^4 erg 


4.7 X 10*V3.5 
4.7 X 1041/3.5 
4.7 X 1041/3.5 
1.5 X 1055/3.5 


3.2 
3.2 
3.2 

3.2, 7.1 


Evidence for pair production 


Bulge 511 kcV Une emission 600 N+ ^ 1.5 x lO^^ 

(production by ee) 
(production by 67) 
(production by 77) 


3.7 X 1042/325 
8.1 X 1041/335 
1.1 X 1042/345 


5.4.1 
5.4.2 
5.4.3 



* X-ray luminosity (Lx), kinetic luminosity Lkin, energy stored in the hot gas (L^iiot), or positron production rate {N+) 
required to explain the observations. 

^Thc dependence on parameters other than the accretion rate is omitted, where /3.5 = /6/103-5 and ft, = A/past / ALprcscnt- 
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Fig. 1. Schematic diagrams for various phenomena on various scales around the Galactic center discussed in this paper. The 
upper- left diagram is the all sky map showing the scales of the 511 keV annihilation line emission and the outflow making the North 
Polar Spur (NPS). The upper-right diagram is for the region within ~ 1° of the Galactic center, showing the scale of the outflow 
observed as the Galactic center lobe (GGL) and the expanding molecular ring (EMR). Famous objects found in the radio image of 
this region are shown by greyscale, and the X-ray reflection nebulae are indicated by the black regions. The lower diagram is for 
the innermost region, showing the interaction of Sgr A East, West, and Sgr A* surrounded by the ionized halo. The gravitational 
radius (dashed circle) is a radius within which the SMBH gravity is dominant compared with the surrounding stars. 



